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Numerical calculations of the thermospheric tidal winds and temperatures were produced by use 
of the NCAR Thermospheric General Circulation Model (TGCM). The present calculations were 
restricted to solar minimum equinox conditions. The effects of viscosity, conductivity, diffusion, 
ion drag, winds, and temperature gradients were included. The semidiurnal propagating waves ex- 
cited by heating in the lower atmosphere were modelled by use of the classical tidal perturbations 
as lower boundary conditions. The TGCM was tuned by adjusting the tidal forcing term until 
calculated semidiurnal winds and temperatures matched incoherent scatter observations. The tidal 
TGCM results are consistent with previous theoretical work and successfully reproduce high 
altitude temperature and velocity data, but give somewhat lower magnitudes for velocities and 
temperatures near 160 km than are seen observationally. Future plans call for similar studies being 
done for solar maximum equinox conditions, and then the solstice cases. These studies would 
benefit enormously from additional data, particularly if higher order modes are to be accurately 
included. 


63 




Figure 1. Classical Hough functions for the semidiurnal tides. Note increasing latitude 

structure for higher order modes. 
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Figure 2. Classical solutions for any 
tidal variable as a function of altitude, 
latitude, and time can be written as 
shown in the top line of the figure. 
The geopotential $ at 97 km can be 
written as shown. 
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1. SPECIFY THE GEOPHYSICAL CONDITIONS; 

2. START WITH EQUILIBRIUM SOLUTIONS OF THE TGCM FOR 
SOLAR HEATING ALONE FOR THE APPROPRIATE CONDITIONS; 

3. GUESS THE AMPLITUDE AND PHASE FOR THE 2,2 AND 2,4 
GEOPOTENTIAL AT 97 KM, AND CALCULATE THE TOTAL TIDAL 
GEOPOTENTIAL AT 97 KM; 

4. CALCULATE U, V, T AT 97 KM FROM THE GEOPOTENTIAL AT 97 KM 
USING FORMULAE FROM CHAPMAN AND LINDZEN AND USE THESE AS 
LOWER BOUNDARY CONDITIONS; 

5. INTEGRATE THE MODEL EQUATIONS UNTIL EQUILIBRIUM IS 

achieved. This typically takes lOh min Cray cpu; 

6. FOURIER DECOMPOSE THE U, V, AND T FIELDS INTO ZONAL 
WAVENUMBERS 1 AND 2; 

7. COMPARE THE TGCM U, V, AND T AMPLITUDES AND PHASES AT 
18° AND 45° WITH OBSERVATIONS FROM ARECIBO, MILLSTONE 
HILL, AND SAINT SANTIN; 

8. MODIFY THE INITIAL GUESSES FOR THE 2,2 AND 2,4 AMPLITUDES 
AND PHASES OF THE GEOPOTENTIAL AND REPEAT FROM STEP 4 
UNTIL REASONABLE AGREEMENT IS OBTAINED BETWEEN THE MODEL 
AND OBSERVATIONS. 

Figure 3. Steps in the calculation of the thermospheric winds and temperatures including the 

effect of tidal excitation. 
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Figure 5. The incoherent scatter radar observations to which the TGCM winds and 

temperatures were compared. 
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BEST FIT TO OBSERVATIONS FOR 

2, 2 AND 2, 4 GEOPOTENTIAL AT 97 KM: 

A 2 2 = 4 x 10 4 m 

A 2 4 = 3 x 10 4 m 

<f> 2 2 ~ “0-5 hr 

</> 2 4 = -2.5 hr 


Figure 6. The amplitudes and phases for the 
2, 2 and 2, 4 geopotential at 97 km which 
best fit the observations shown in Figure 5. 
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Figure 8. The TGCM north-south velocity at 18° and 45°. 
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AMPLITUDE (°K) HOUR OF MAX 

Figure 9. The temperature observations at 18° and 45°. 
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Figure 10. The TGCM temperatures at 18° and 45°. 
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OBSERVATIONS 



Figure 1 1 . The east-west velocity observations at 45 
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Figure 12. The TGCM east-west velocities at 18° and 45°. 
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Figure 14. The TGCM diurnal north-south velocity results at 45 °. 
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